Microscopic origin of critical current fluctuations in large, small and ultra-small area 

Josephson junctions. 
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We analyze data on the critical current and normal state resistance noise in Josephson junctions 
and argue that the noise in the critical current is due to a mechanism that is absent in the 
normal state. We estimate the noise produced by conventional Two Level Systems (TLSs) in the 
insulating barrier and find that it agrees both in magnitude and in temperature dependence with 
the resistance fluctuations in the normal state but it is not sufficient to explain the critical current 
noise observed in large superconducting contacts. We propose a novel microscopic mechanism for 
the noise in the superconducting state in which the noise is due to electron tunneling between 
weak Kondo states at subgap energies. We argue that the noise produced by this mechanism gives 
temperature, area dependence and intensity that agree with the data. 
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Introduction.— The microscopic mechanism at the ori- 
gin of critical current fluctuations in Josephson junc- 
tions is a long standing open problem. A phenomeno- 
logical characterization of critical current noise began 
with the experiments in large Niobium junctions in the 
late 80's in these experiments the critical current 
fluctuations were inferred from the fluctuations in the 
current of resistively shunted junctions biased by a low 
voltage that has no effect on the superconductivity of 
the contacts themselves. It was found that the noise 
power spectra display 1// behavior for low frequencies 
/ < lOHz with r^— temperature dependent intensity 
and that it also scales with the inverse area of the junc- 
tion with coefficient that depends, albeit weakly, on the 
superconducting or insulating barrier material [3| . 

It is natural to attribute the critical current fluctu- 
ations to charges that move between different localized 
states in the junction barrier (3]. However, a more de- 
tailed comparison with the experiments reveals serious 
problems with this picture. In this model, each fluctu- 
ator is similar to a glassy TLS 0], so one expects that 
the distance between the levels and the tunneling barrier 
have a broad distribution on the atomic energy scale. 
Quantum tunneling or thermal activation leads to the 
charge motion between two states which changes the bar- 
rier height and thus leads to a noisy Random Telegraph 
Signal (RTS) in the current through the junction [6|]. A 
superposition of RTSs with a broad distribution of tun- 
nehng barriers leads to the 1// noise power spectra 0. 
The problem with this model is that TLSs and similar 
objects have a constant density of states at low energies 
and this would give a linear T— dependence of the noise 
power spectrum in contrast with the data. Alternatively, 
assuming that behavior comes from the linear den- 
sity of states and extrapolating this density of states to 
atomic scales one gets unphysically large values for the 



density of these switches thus indicating the presence of 
some low energy scale in the problem. 

Recently, the interest in critical current fluctuations 
was renewed because of their importance to the super- 
conducting qubit dcphasing and a new puzzle was added 
to the picture. The new experiments Q studied the fluc- 
tuations in the small area {A ^ O.lfim^) Aluminium junc- 
tions, similar to the ones that arc currently used for qubit 
implementations in the flux qubit Q , the phase qubit [lo| 
and the quantronium [llj . The critical current noise was 
inferred from fluctuations in the normal state resistance 
of the junction. It has been observed that the noise power 
spectrum displays a 1// behavior at low frequency but 
its temperature dependence is only linear and the inten- 
sity of the noise is two order of magnitude lower than the 
value reported for larger superconducting contacts. 

In this Letter we argue that the problems with the mi- 
croscopic model noted above as well as the inconsitency 
between the old and the new data are all removed if the 
main source of the critical current noise is electron trap- 
ping in shallow subgap states that might be formed at 
the superconductor- insulator boundary The elec- 

tron tunneling between such traps contains two Fermi 
factors leading to the T^/Tq behavior with the energy 
scale To ~ Tc, thereby eliminating the difficulties with 
the conventional model noted above. This model is fur- 
ther supported by the recent experiments performed in 
Single Electron Transistors that show dependence of 
the low 1// charge noise power spectra in very small 
superconducting contacts [13[. We shall argue that this 
mechanism dominates the critical current noise in the su- 
perconductive regime but disappears in the normal state 
where these states become Kondo resonances leaving only 
a weaker conventional TLSs mechanism active that, at 
high temperatures, produces linear T— dependence of the 
noise power in agreement with the data Q . An important 
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feature of the new mechanism is that the large number 
of subgap states is partially compensated by their small 
weight (a vestige of their Kondo resonance origin) [l^ . 
Physically, it means that these states generate a feature- 
less 1// noise even in the smallest contacts. In contrast, 
in normal leads, where the noise is due to TLSs in the 
insulating barrier, our estimates show that the number 
of active fluctuators, especially in ultrasmall contacts, is 
low, in agreement with the data 

Below we discuss the details of the data and estimates 
of the noise produced in new and conventional micro- 
scopic models. We begin with the data. 

Data.— The directly measured low frequency noise 
(/ < lOHz) of the critical current turns out to be almost 
universal. The original experiments measured it in large 
contacts (area. A, varying between 10 — 100/zm^) made 
of Nb-NbOx-PbIn P, Hi- In these experiments the junc- 
tions were biased by a very low voltage V Ki 1 — BfiV that 
does not affect the superconductivity of the contacts (i.e. 
y<C 2A/e) [l^. In these conditions, the noise power 
spectra are described by: 
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in a wide temperature range 0.09K <T< 4.2K. Here /o 
is the critical current, Aq and Tq are the area and tem- 
perature scale which are conventionally set to Aq = 1/^m^ 
and To = 4.2K. Remarkably, the dimensionless propor- 
tionality coefficient, 7, is not sensitive to the details of 
the junction preparation and it is roughly (within a factor 
3) "universal ": 7 « 1.44 • 10'^° The noise observed 
in these experiments was usually featureless, evidently 
coming from many fluctuators, but in rare cases and at 
relatively high temperatures, one observes also a switch- 
ing process between two well defined current values that 
disappears when temperature decreases below T < IK. 
Notice that the quadratic growth of the noise with tem- 
perature that starts at T = O.IK implies that the num- 
ber of fluctuators at atomic energies would be at least 
lO^'' larger than their number at low temperatures and 
thus would exceed the number of atoms in 10/xm^ contact 
{N < 10^), providing the evidence of additional energy 
scale. 

In a different set of experiments Q, the critical cur- 
rent of small junctions was not directly measured but ex- 
tracted from their normal state resistance fluctuations. 
To perform measurements at very low temperatures, the 
Al-AlOaj-Al junctions of ^ ^ 0.1/xm^ area were subjected 
to the magnetic field above O.IT that suppressed the 
superconductivity in the contacts. It was found that 
for 0.3K <T< 5K the noise power spectrum is well de- 
scribed by a linear T-dependence: 
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with the intensity that is two order of magnitudes lower 
than that predicted by Eq. ([J). Moreover, for tempera- 
tures below 0.8K, it has been observed that the low fre- 
quency noise was due to few individual strong fluctuators. 
These observations are in a good agreement with the con- 
ventional TLSs picture of the noise origin. Indeed, in a 
typical structural glass one expects a constant density of 
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eV that leads to the linear tem- 



perature dependence of the noise and to the total num- 
ber of the thermally excited TLSs, A/'TLSs(r = IK) ~ 10, 
in contact. These data are also in a qualitative 

agreement with the older measurements of ultrasmall 
(A < 0.05/im^) Nb-Nb205-Pbln junctions performed in 



the temperature range IK <T < 300K [lj|. In these 
measurements, the contacts cither were kept normal due 
to a high temperature or they were very far from equilib- 
rium due to a high voltage bias V > lOmV ^ 2A/e; the 
observed low frequency 1// noise was dominated by few 
fluctuators even at T < lOOK. 

The detailed quantitative comparison between differ- 
ent data and microscopic estimates is made difficult by 
the possibility of the inhomogeneous current distribution 
through the insulating barrier. For instance, the newer 
data [SI] report the relative variation of the resistivity due 
to one fluctuator corresponding to the change in the ef- 
fective area of the junction by SA^ff = Inm^ while for 
the older data [3l this value is SAejf lOnm^. This 
difference explains why the noise in new experiments got 
dominated by a few fluctuators only below IK in con- 
trast to where single fluctuators were resolved below 
lOOK. Indeed, it is difficult to imagine microscopic pro- 
cesses that affect large areas ~ lOnm^; a more plausible 
explanation is that the conductance is dominated by rel- 
ative small areas or channels in this system [l6| which 
decreases further the number of fluctuators (Figure 1). 
This conjecture is supported by the recent analysis [17| 
that found that the conductive area is At ~ 10~^ of the 
total. Most likely, however, this value depends strongly 
on the material preparation. 
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FIG. 1: Sketches of the microscopic mechanisms of decoher- 
ence that might be responsible for critical current noise. 

Microscopic mechanisms of decoherence.— All mech- 
anisms of the noise assume that it comes from the 
fluctuations in the state of the insulator that affect 
the tunneling barrier and thus the critical current. If 
these fluctuations are not affected by the state of the 
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metal (superconductor or normal) they would equally af- 
fect the critical current Iq and the normal state resis- 
tance Rn according to the Ambegaokar-Baratoff formula: 
lo = 2^A(T)tanh (^^^^ [Hi. One can argue that 

these fluctuations might be due to fluctutating trapped 
electrons or randomly moving TLSs inside the barrier. 

Low energy electrons trapped in the vicinity of the nor- 
mal metal can tunnel into the metal; this process converts 
the sharp state into a short-living resonance. In super- 
conductor, however, the low energy states disappear and 
the life time of these states increases. Notice that a large 
potential difference between the contacts provides elec- 
tron states below the gap and thus makes the system 
effectively normal. In a superconductor a localized low 
energy electron can change its state only by tunneling 
into another subgap state via Andreev process [l^. If 
the density of states for these subgap states is constant, 
the number of thermally active electrons and their possi- 
ble final states both scale with temperature leading to 
behavior of the noise in agreement with the data. How- 
ever, similarly to TLSs, the density of states for such 
traps is too little to account for the number of fluctua- 
tors observed experimentally. 

We believe that the physical mechanism responsible for 
a dramatic increase in the density of states for the local- 
ized low energy electrons is the formation of Kondo-like 
resonances due to a large Coulomb repulsion between the 
electrons in the same trap [12|. In this picture, the char- 
acteristic energy scale for these resonances is given by 
the Kondo temperature Tk which depends exponentially 
on the bare level width T and the bare level position eo: 
Tk oc exp(— 7reo/2r). The natural assumption that F is 
distributed in a broad range leads to P{Tk) oc \/Tk dis- 
tribution of Kondo temperatures in the normal state. In 
the superconductor the resonances having « 0.3A be- 
come localized low energy levels with a constant surface 
density of states z^(e) = p2D /(^Tk) where p2D is the bare 
surface density of traps {p2D~ 10^^/cm-^) and *p w 1. 
The increase in the density of states of these objects is 
partially offset by the small weight of the Kondo reso- 
nance w ~ Ta'/eq that translates into the small weight of 
the formed localized state. 

The noise at frequency / is due to the electron tun- 
neling between two traps with rate r^^ ^ f. The exact 
calculation of the rate for two specific traps is a compli- 
cated problem but it is sufficient for our purposes to note 
that the latter is an exponential function of the bare pa- 
rameters and thus it is characterized by dPlr) = dr/r 
distribution function. Indeed, in superconductors the 
tunneling amplitude falls exponentially with distance: 
t oc e~^^^ where f is the coherence length of the super- 
conductor. Moreover, the tunneling between Kondo res- 
onances is additionally suppressed by the small factor 
w and finally, the electron tunneling between the states 
with different energies is also accompanied by the ther- 



mal phonon emission or absorption but the rate of these 
processes is much faster than the frequencies of interest 
(/ < IkHz) so the main dependence comes from the 
exponential factors mentioned above. 

In order to estimate the critical current fluctuations in- 
duced by these tunneling processes we need to know the 
effect induced by a single electron trapped in a deep state. 
We shall describe it by the change in the effective area of 

the junction SA^ff (i.e. SIq = -^SAeff). For the quan- 

titative estimates we shall assume that SAef / ~ ^Aej / . 
The tunneling processes allow tunneling between the 
traps within energy T from the Fermi surface located at 
distances r < ^ from each other. Combining all factors 
we find that the noise spectrum generated by indepen- 
dent relaxational processes between the traps reads: 
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where c ~ 0(1). This noise displays dependence and 
inverse proportionality to the area of the contacts in 
agreement with the data observed in large contacts. To 
estimate the intensity of the noise we assume that the dis- 
tance between electron traps is roughly of the same order 
as the distance between TLSs in the bulk or (between 
typical surface defects) leading to p2D ~ lO^^cm"^. 
Using the Nb parameters: ANb = 16.33K, ^Nb ~ 40nm, 
SAeff ~ 5nm^, w = 10~^, we estimate the dimensionless 
parameter 7 controlling the noise intensity at T = 4.2K, 
obtaining 7 « 5 • 10^ which is in excellent agreement 
with the measured value (given the crude nature of the 
estimates). 

Alternatively, critical current noise might be due to 
fluctuating TLSs in the insulating barrier [5| . Most likely, 
TLSs correspond to atoms that can tunnel between two 
positions in the amorphous material. Two different states 
can be distinguished by the different values of the dipole 
moment. Thus, TLSs deep inside the insulator inter- 
act only weakly with electrons in the metal and are not 
affected by the superconductivity. Each quantum TLS 
is described by the Hamiltonian Hq = eaz + tax , where 
e is the energy difference between the two minima and 
t = hojQe'^ is the tunneling amplitude between them. 
Here luq is of the order of the frequency vibrations of the 
particle in the potential wells and S* ^ 1 is the tunneling 
action. In bulk amorphous materials e and S are dis- 
tributed in a broad range leading to the distribution func- 
tion P{e,t) ~ v/t. Generally, each TLS has states with 
energies E = ±y/e'^ + that are occupied with ther- 
mal probability factors. At not-too-low temperatures 
(T > 0.2K) the TLS relaxation process is dominated by 
phonon emission or absorption that leads to the relax- 
ation rate of the occupancy number r^^^ = aEt^ coth [3E 
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where the coefficient a depends on the details of the 
TLS-phonon interaction [20|. The dominant contribu- 
tion to low frequency noise comes from thermally ex- 
cited TLSs with low tunneling amplitudes t ^ e ^ T. 
These TLSs are characterized by the distribution func- 
tion P{E,t) ~ vjlT. Each individual TLS contributes: 

s'^-f^ = — 1 I Tii ■> to the current-current correlator. 
By assuming as before that each TLS affects the criti- 
cal current by SIq ~ {SA/A)Io and by averaging over the 
TLSs distribution function wc find that the total noise 
power spectrum originating from the uniform insulating 
barrier of thickness d reads: 

SiJu;)::^ K^SA^iyd-. (4) 

°^ ' A to ^ ^ 

As expected this noise is proportional to T and 1/A. Non 
uniform distribution of the current discussed above de- 
creases the number of the effective TLSs because each 
TLS affects only the current that flows in a close vicinity 
of a moving atom. We describe this suppression by the 
dimensionless coefficient k; very crudely we can estimate 
k^6Aq/5A where JAq = O.lnm^. For small Al con- 
tacts, A = Ifim"^, at T = 4.2K we get SiJI^ « 3 • 10''^^, 
which is in agreement with the direct data on the nor- 
mal state resistance fluctuations 0] . Notice that a rel- 
ative small value of the effective area SA for these con- 
tacts implies that the current is relatively homogeneous; 
one expects smaller values for the contacts studied ear- 
lier [14| for which similar estimates gives a larger result 
Sig/Ig ~ 3 • 10~^^ which is however still smaller that the 
noise observed directly in the superconducting state. 

As a final remark, let us notice that the electron trap 
model for the noise in the superconducting state leads 
to an interesting prediction. Namely, the process of tun- 
neling between these traps should also contribute to the 
high frequency dissipation leading to ohmic behaviour 
of the spectrum: 5*70(0;) = 2e^Gui. The proportionality 
constant, G, in this formula has a meaning of the effective 
conductance of the contact due to the trapped electrons, 

i.e. G « 7^ ( 2e'^R To ) ' could be in principle ex- 

tracted from measurements similar to the ones performed 
by Astafiev et al. [l^l for the charge noise. Here, how- 
ever, such experiments are made difficult by the noise 
coming from the shunt resistors and the residual resis- 
tance of the contact itself due to charge fluctuators. 

Conclusions. — We identify a new microscopic mech- 
anism responsible for the noise in the superconducting 
contacts. This mechanism contributes to the low fre- 
quency noise only in the superconducting state and it is 
due to the appearance of subgap states localized near the 
superconductor-insulator boundary. The resulting noise 
has homogeneous spectrum, dependence and it is in- 
versely proportional to the surface area A of the junction 
in agreement with the experiments. The estimates of its 
intensity are in a good agreement with experimental val- 



ues. This mechanism disappears when the superconduc- 
tivity of the contacts is suppressed. In the normal state 
the noise is generated by the thermally excited TLSs. In 
this case, the noise power spectra of the critical current 
fluctuations show a linear T dependence, inverse propor- 
tionality to the area, with intensity of the noise lower 
than the one generated by the superconducting state in 
agreement with experiments. The mechanism that dom- 
inates the critical current noise should also dominate the 
charge noise. It would be important to confirm it directly 
by measuring the charge and critical current/resistance 
noise on the same (or at least similar) samples in both 
superconducting and normal states. 
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